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We explore the dynamics of assisted quintessence, where more than one scalar field is present 
with the same potential. For potentials with tracking solutions, the fields naturally approach the 
same values — in the context of inflation this leads to the assisted inflation phenomenon where 
several fields can cooperate to drive a period of inflation though none is able to individually. For 
exponential potentials, we study the fixed points and their stability confirming results already in 
the literature, and then carry out a numerical analysis to show how assisted quintessence is realized. 
For inverse power-law potentials, we find by contrast that there is no assisted behaviour — indeed 
those are the unique (monotonic) potentials where several fields together behave just as a single 
field in the same potential. More generally, we provide an algorithm for generating a single-field 
potential giving equivalent dynamics to multi-field assisted quintessence. 
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I. INTRODUCTION 

An attractive hypothesis for the fundamental nature 
of dark energy is quintessence — a scalar field evolving 
in a non-zero potential energy 0. Such modelling has 
proven highly successful in implementing inflation mod- 
els in the early Universe. It has however met with much 
less success in the present Universe, primarily due to the 
difficulty of obtaining equations of state close enough to 
the cosmological constant value w = — 1 to satisfy obser- 
vational bounds 0- One strategy is to ensure that the 
quintessence energy density starts at such a low value 
that the field only begins evolving close to the present, 
but such tuning is hardly more satisfactory than a cosmo- 
logical constant. The alternative is for the field to evolve 
significantly, ideally exploiting a 'tracker' behaviour ren- 
dering the late-time evolution almost independent of ini- 
tial conditions. Unfortunately however this is viable only 
for very particular potentials: amongst monotonic poten- 
tials, exponentials do not give acceleration and power- 
laws can give negative enough w only for exponents well 
below 2. Among the potentials regarded as giving satis- 
factory phenomenology, many in fact feature a minimum 
tuned to match the observed value of the cosmological 
constant. 

In this paper we investigate whether this situation 
might be alleviated by allowing the quintessence to arise 
from several fields, which we assume to have the same 
potential energies. While our interest is primarily phe- 
nomenological, we note that such situations may arise 
from higher-dimensional theories; in fact there are many 
dynamical modulus fields in string theory corresponding 
to the size of compactified dimensions. In the context 
of early Universe inflation, such a set of fields has been 
shown to give the phenomenon of assisted inflation 0, 
whereby they may collectively drive inflation even if each 
individual field has too steep a potential to do so on its 
own, i.e. yielding an effective equation of state closer to 
W = — 1. It is therefore worth considering the possi- 



bility of assisted quintessence behaviour, in order to see 
whether it may be better able to match observations. As- 
sisted quintessence is particularly attractive in the con- 
text of tracking models, as each field will separately con- 
verge onto the tracking solution making it entirely natu- 
ral that they all play a dynamical role. 

In this paper we will investigate some aspects of as- 
sisted quintessence, mostly restricting ourselves to the 
simplest case where each field has the same potential and 
always assuming there are no interactions between fields. 
The most closely-related paper is that of Blais and Po- 
larski who analyzed several multi-field quintessence 
models both analytically and numerically, though with 
a different focus directed mainly at attempting to re- 
alize models where the acceleration is a transient phe- 
nomenon. We will primarily study the exponential and 
inverse power-law cases in detail, in the former case con- 
sidering different exponents for the two potentials as al- 
ready extensively analyzed by Coley and van den Hoogen 
0. In this paper we will also provide an algorithm for 
relating assisted quintessence to an equivalent single-field 
model under more general circumstances. 

Our scenario is distinct from two types of scenario al- 
ready extensively investigated in the literature. One is 
the double exponential potential models of Refs. 
where there was only a single field (Ref. |i| briefly 
mentioned a multiple-field case but not with uncoupled 
fields). Another is the two-field models of dark energy 
which have received some attention recently as a way of 
crossing the cosmological constant boundary (w = — 1) to 
give rise to a phantom behaviour at late times 0, 0, [h| . 
Unlike our case, at least one of those fields must have a 
negative kinetic energy. 



II. ASSISTED QUINTESSENCE 

If one accepts the possibility of multiple fields, par- 
ticularly sharing the same form of potential V(4>i), then 
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the idea of assisted quintessence emerges very naturally 
provided the potentials have tracker solutions for at least 
some values of the fields. Tracker solutions arise when the 
fields are initially sub-dominant as compared to a perfect 
fluid, the field contribution to the Friedmann equation 
then being neglected to give equations of the form 
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where H is the Hubble rate, mpi is the Planck mass and a 
dot denotes the derivative with respect to a cosmic time 
t. Here the fluid density pt might for instance be broken 
up into matter and radiation components p nl and p T . In 
this set-up, the scalar fields are completely unaware of 
each other's existence (their normal channel of commu- 
nication being via the Friedmann equation), and hence 
separately evolve onto the tracker solution in response to 
the fluid. At sufficiently-late times one would therefore 
have all the <pi equal to each other, and hence of equal 
potential energy. 

For an individual field <f> with potential V(<fi), whether 
or not there is tracking behaviour can be determined from 
the value of the function 

VV" 

where a prime represents the derivative in terms of ip. 
Solutions converge to a tracker provided that it satisfies 
r > 1 — (2 — 7)/ (4 + 27) where the equation of state is 
Pf = (7 — l)p{, the interesting case however being T > 1 
which is required for the field energy density to grow 
relative to the fluid allowing eventual domination 11]. 

The convergence of different fields to the same track- 
ing solution does not in itself amount to assisted 
quintessence, as the fields are not generating any grav- 
itational effect on the background evolution. The con- 
vergence does however set up the initial conditions for 
such a behaviour. What is mainly of interest is what 
happens once the energy density of the fields, all evolv- 
ing together, is no longer subdominant, and that is the 
situation addressed in the rest of this paper. We will 
therefore be considering the full Friedmann equation 
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where p ( f >i = Vi((f>i) + (f>f/2 and we assume spatial flatness 
throughout. 

III. EXPONENTIAL POTENTIALS 



where k 2 = 87r/m P1 . For generality, we will allow the 
potentials to have different slopes. 



A. Assisted quintessence solutions 

For the case where no matter is present, this system is 
exactly the original assisted inflation scenario of Liddle et 
al. @, where the multiple fields evolve to give dynamics 
matching a single-field model with 
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For a single-field potential V(4>) — Voe^ XK ^ the scale fac- 
tor evolves as a k F, where p = 2/ X 2 . Then the multi- 
field case Eq. (0 corresponds to an effective power-law 
index given by p e g = 2/ X 2 S = ^2pi • The expansion 
rate is therefore more rapid the more fields there are. 
A particularly comprehensive analysis of multiple fields 
in exponential potentials has been given by Collinucci et 

al. H3. 

As it happens, this assisted behaviour continues to be 
valid in the presence of a perfect fluid with equation of 
state p = (7 — lYp. This was first noted by Coley and 
van den Hoogen jfj, who also allowed for the possibility 
of spatial curvature. This is because the method used to 
relate multi-field dynamics to an equivalent single-field 
dynamics in Ref. 3] is a property of the scalar field sector 
alone, being valid in the presence of any other matter 
sources. 



B. Critical points and stability 

The critical points for this system were classified by 
Coley and van den Hoogen [B| , and we will reiterate their 
results only briefly before embarking on some numerical 
analysis of the evolution towards them. In Ref. spa- 
tial curvature was included as a degree of freedom and 
the fluid assumed to have 7 > 1, whereas we assume 
spatial flatness and 7 in the wider range to 2. 

To study the critical point structure and stability of 
the system, it is convenient to introduce the following 
dimensionless quantities [13J 



We consider two fields <pi and 4>2 each with a separate K X i K\fVl 

exponential potential x i — ^/qjj' ^ i ~ yfojj ' * = 1) 2 • (7) 

V((t>i,<f> 2 ) = Ae-^ 1 + Be-^ 2 (5) 

= Vi((/>i) + V2(<fo) i Then we obtain 
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TABLE I: The properties of the critical points in the presence of a barotropic fluid with < 7 < 2. There are seven discrete 
points, while Case 2 corresponds to a circle of points parametrized by an angle 6. Case 2 is neutrally stable along the circle, 
but always unstable in at least one other direction. The last two columns show the effective energy density and equation of 
state of the two scalar fields combined. 
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-3xi + \J^Vi + [2a;? + 2a; 2 , + 7(1 - a;? - y\ - x\ - y\)\ 
l~3 3 

-AiW - Xjj/j + -yi [2x\ + 2x\ + 7(1 - x x - y x - x 2 - y 2 )] 
-| [2a; 2 + 2a; 2 + 7 (1 - a; 2 - y 2 - x\ - y 2 )] , 



i = 1,2, 
i = 1,2. 
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where = In a, together with the constraint 

(11) 

and the equation 
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We define the density parameters fi^ 
of state Wfc for scalar fields, as 



i = 1,2. 



■vt 



(12) 



We note that the density parameter for matter is given 
by Q m = K 2 pt/3H 2 = 1 - - fl^ from Eq. l(TT|l . 

These equations are in fact valid for any uncoupled 
potentials, with A^ defined by 



A,- = - 
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(13) 



Our case of exponential potentials corresponds to A,; both 
constant. It is straightforward to extend our analysis to 
the case of a dynamically changing A as studied in single- 
field models in Ref. Q. 

The classification of critical points and their stability 
is given m Tabled and agrees with results in Ref. Q. 1 



There are eight types of critical point, seven being dis- 
crete points and one a circular locus of critical points. 
They are readily compared to the five types of critical 
point found for the single-field system in Ref. |l3|. Our 
cases 1, 2 with 9 = and it, 3 and 4 correspond to the 
five points in the single-field case for the field 4>i, with 
4>2 playing no role. Our cases 1, 2 with 9 = n/2 and 
3tt/2, 5 and 6 are the same solutions for the 4>2 field. Fi- 
nally, cases 7 and 8 are critical points in which both fields 
play a role, and which have no direct analogue with the 
single-field case. 

In the single-field case, the stable late-time attractor 
is either scalar-field dominated (case 3 or 5) or a scaling 
solutionjcase 4 or 6) depending on the relative values of A 
and 7 [l3| . Once a second field is added, the new degrees 
of freedom always render those solutions unstable. The 
late-time attractors instead become either the assisted 
scalar-field dominated solution case 7 (for A 2 ff < 3 7 ) or 
the assisted scaling solution case 8 (for A 2 ff > 37). 

To compare our results with Ref. Q , we note that un- 
der their assumptions the curvature always dominates the 



We note that a full analysis has also been carried out for the case 



where one of the fields is of phantom type (opposite sign of the 
kinetic energy) [jj], where the late-time behaviour is domination 
by the phantom field. 
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FIG. 1: A numerical integration for the case Ai=2, A2 = 1.8 
and 7=1, with initial conditions chosen so that <j>\ domi- 
nates. Field 1 is in solid and field 2 dashed. The solution 
initially approaches the single-field scaling solution (case 4), 
but this becomes unstable once the second field becomes im- 
portant, the late-time attractor being the multi-scalar-field 
dominated accelerating solution (case 7, A e g = 1.34). The 
dotted line shows — q = a/aH 2 , where q is the deceleration 
parameter. 



fluid at late times, and behaves like a 7 = 2/3 fluid. Ac- 
cordingly, they always find solutions with non-zero fluid 
density to be unstable to eventual curvature domination, 
whereas our assumption of spatial flatness renders them 
stable. Otherwise, the classification and stability shown 
in our Table matches their Tables I and III. Addition- 
ally they correctly describe the spatially-flat case in their 
text. 



C. Multi-field phenomenology 

For the exponential potential, the assisted quintessence 
phenomenon does indeed exist, with the extra fields be- 
ing equivalent to a flatter single-field potential. Unfortu- 
nately this result does not seem particularly useful phc- 
nomenologically, as the scaling solutions do not give ac- 
celeration as required by observations, while the scalar- 
field dominated solutions would long ago have made the 
matter density negligible. 

There is however one scenario that might be of inter- 
est, which is to imagine there are a large number of ex- 
ponential potentials with different initial conditions. As 
the Universe evolves, more and more fields would join 
the assisted quintessence attractor, reducing A g. Even- 
tually, this could switch the attractor from the scaling 
regime A^ ff > 37 into the regime of late-time scalar field 
dominance Ag ff < 37 pi- 
Figure^ shows a two-field example of such evolution, 
with the single-field scaling solution becoming unstable 
as the second field becomes important, switching the 
evolution into late-time scalar field domination. From 
Table I the first scaling regime corresponds to case 4 
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FIG. 2: A comparison of two-field evolution for Ai = 3 and 
A2 = 1 (top) and one-field evolution for A = 1 (bottom) . The 
density parameters Q. ( f, i and Q, m are shown in solid lines and 
the equation of state parameters Wfc and w in dashed lines. 
The region around 6 < N < 11 corresponds to the scaling 
regime in which the energy density of the field <j>i decreases 
proportionally to that of the fluid. The present epoch is de- 
fined in each simulation as the time when the matter density 
drops to 0.3 of critical. 



[x\ = yi = y/3/8, X2 = 1/2 = 0) with no acceleration, 
whereas the final stable attractor is case 7 {x\ — 0.365, 
yi = 0.560, x 2 — 0.406, 1/2 = 0.623) with acceleration. 
We checked that the values of the fixed points agree very 
well with numerical results. 

However even if the initial conditions were fine-tuned 
to bring this transition into the recent past, it is hard to 
see how the equation of state w could reach a sufficiently- 
negative value to be observationally viable, current ob- 
servations indicating roughly w E7-I < -0.8 @. In 
fact, if we compare the two-field scenarios with a single- 
field scenario in which an accelerated expansion occurs at 
late times (A < \/2), we find that typically the two-field 
models give a larger present value of w (here w is the 
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equation of state of scalar fields including the contribu- 
tion of both 4>i and </> 2 ). Figure |3 compares a two-field 
model with Ai = 3 and A2 = 1 against a single-field 
model with A = 1, in both cases with the field with slope 
1 starting with a low value and coming to dominate only 
at the present epoch (the 'thawing' regime as according 
to Ref. We see that the equation of state of the 

thawing field alone is indeed closer to —1 at the present 
(identified as when the matter density is 0.3 of critical) 
than in the single-field case, but the combined w of the 
two fields is larger. 



IV. INVERSE POWER-LAW POTENTIALS 

We now consider the inverse-power law case V((f>) = 
Vo</> - ' 3 , where we will take the same exponent and nor- 
malization for each field. Here one might hope that 
the assisted phenomenon would give rise to an effec- 
tive P e g which is smaller than the individual /3, so that 
observationally-viable models can be achieved in steeper 
potentials than the single-field case. Unfortunately that 
turns out not to be true, as we now see. 

These potentials are favoured because they exhibit 
tracker solutions, and it is therefore legitimate to suppose 
that after some early time we can take <px = (f>2 = ■ ■ . = (f>. 
The Friedmann and scalar wave equations become 



behaviour so that we can take 



giving 



H 2 = 



8tt 
3mpj 



Pi + nVo4> 



= -3H + pVatf)- - 1 



(14) 
(15) 



where n is the number of fields. This is not yet equiva- 
lent to a single-field model, but can be made so (loosely 
following the method of Ref. [3j]) by the redefinitions 



X 



W = VnP nV 



(16) 



which then gives the equations of a single-field model 
with potential W(x) = Wox -13 - 

This rescaling indicates that the assisted behaviour 
renormalizes the amplitude of the potential in this case, 
but does not renormalize its exponent. In such models 
the amplitude has to be adjusted in order to give the 
present-day value of the matter density f2 m , and having 
done that the multi-field system then has identical dy- 
namics to the single-field model, in particular predicting 
the same present-day value of the equation of state w. 



V. ASSISTED QUINTESSENCE DYNAMICS 

We end with a more general construction for analyzing 
assisted quintessence, extending the analysis of the pre- 
vious section to the case where the potential is arbitrary, 
but the same for all fields. Again we assume tracking 
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(17) 
(18) 



To find an equivalent single-field system, we first note 
that the kinetic term in the Friedmann equation forces 
the correspondence 



(19) 



The equations can then be transformed into single-field 
form with potential 



W( X ) = nV( X /Vn~) , 



(20) 



which obviously has the desired effect in the Friedmann 
equation, but which also renders the fluid equation into 
single-field form. This formula therefore represents an 
algorithm for finding a single-field potential W which will 
generate the same evolution as multiple fields evolving 
together in the potential V. 

For simplicity, we henceforth consider the two-field 
case, though the generalization is straightforward. The 
correspondence then is 



X = \/20 ; W(x) = 2V( X /V2) 



(21) 



This allows us to ask what condition would have to be 
satisfied in order to have no assisted behaviour, i.e. for 
W and V have the same functional form apart from an 
overall constant. This happens for potentials obeying the 
condition 



V(x) = 2CV( X /V2) , 



(22) 



for all x, where C is a constant. The general solution to 
this equation is 



V( X ) <x x fix) , 



(23) 



where f(x) is any function periodic in lax with period 
ln\/2 (i.e. a Fourier series with this periodicity). While 
this can be any of an infinite class of potentials, having 
such a periodicity is clearly artificial. The only inter- 
esting case therefore is f(x) equals a constant, giving 
the power-law potential. This proves that the inverse 
power-law potentials are the unique monotonic potentials 
which do not exhibit assisted quintessence behaviour. In 
all other cases, the equivalent single-field potential has a 
different functional form. In the exponential case this al- 
gorithm correctly reproduces Eq. © for the case Ai = A 2 . 

As a final point, we note that under the correspondence 
Eq. the tracking parameter T w = WW" /W' 2 is 

equal to Ty = VV" /V' 2 (the primes here being deriva- 
tives wrt the arguments x an d </> respectively). As one 
would expect, the tracking conditions on the multi- field 
model and its single-field dynamical equivalent are the 
same. 
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VI. CONCLUSIONS 

We have studied various aspects of assisted 
quintessence dynamics. Such dynamics arises natu- 
rally if there are several fields with the same potential, 
provided the potential exhibits tracking behaviour for at 
least some stage of its early evolution. 

Our most powerful result is Eq. (j2t)(l . which pro- 
vides a general algorithm for finding a single-field model 
which mimics the dynamics of a multi-field assisted 
quintessence model. Applied to inverse power-law mod- 
els, it shows that they are the unique (monotonic) po- 
tentials for which there is no assisted behaviour, the col- 
lection of fields behaving as a single field in the same 
potential (up to overall normalization). All other po- 
tentials will exhibit assisted behaviour, the exponential 
potential being an explicit example. 

It should be possible to extend our analysis to the case 
of more general dark energy models in which the La- 
grangian includes non-canonical kinematic terms, such 
as the tachyon, k-essence and ghost condensate. For the- 
ories whose Lagrangian p is a function of the field <j> and 



X = (V0) 2 /2, the existence of scaling solutions restricts 
the form of Lagrangian to be p — Xq{Xe x ^), where g is 
an arbitrary function and A is a constant Hi . It would 
be certainly of interest to investigate whether the assisted 
behaviour we found for the canonical scalar field with an 
exponential potential persists in such general dark energy 
models. 

It is interesting that the multi-field system can be an- 
alyzed so simply. Regrettably, however, we have not 
uncovered any scenarios where the assisted quintessence 
phenomenon appears to improve the situation with re- 
gard to the observations. In fact our results show the 
contrary; in potentially the most interesting scenario of 
the inverse power-law it turns out that there is no assisted 
quintessence effect. 
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